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ABSTRACT 

We analyze a complete sample of 7-ray burst afterglows, and find X-ray evidence for high column 
densities of gas around them. The column densities are in the range 10 22 — 10 23 cm -2 , which is right 
around the average column density of Galactic giant molecular clouds. We also estimate the cloud sizes 
to be 10-30 pc, implying masses <; 10 5 M©. This strongly suggests that 7-ray bursts lie within star 
forming regions, and therefore argues against neutron star mergers and for collapses of massive stars as 
their sources. The optical extinctions, however, are 10-100 times smaller than expected from the high 
column densities. This confirms theoretical findings that the early hard radiation from 7-ray bursts and 
their afterglows can destroy the dust in their environment, thus carving a path for the afterglow light out 
of the molecular cloud. Because of the self-created low extinction and location in star-forming regions, 
we expect 7-ray bursts to provide a relatively unbiased sample of high-redshift star formation. Thus 
they may help resolve what is the typical environment of high-redshift star formation. 

Subject headings: gamma rays: bursts — cosmology — extinction — stars: formation 



1. INTRODUCTION 



et al. 2000 ). Despite their unusually high luminosity these 



7-ray bursts emit up to 10 53 ergs in gamma rays and 
X rays in a few seconds, followed by an 'afterglow' of 
X-ray ( Costa ct al. 1997|), optica l (Van Paradijs et al. 
1997) and radio (Frail ct al. 1997) emission that gener- 



ally la sts da ys to months (Van Paradijs, Kouveliotou, & 
Wijers 200C). The arcsecond localizations of 7-ray bursts 
by the detection of these counterparts have made it possi- 
ble to study the environments in which 7-ray bursts arise. 
This has provided a number of indications of their as- 
sociation with massive stars and star formation. First, 
most 7-ray bursts lie within the region of UV emission 
from massive stars in their host galaxy. Second, the en- 
ergies of 7-ray bursts are comparable to those of super- 
novae, further suggesting the deep gravitational collapse 
of a few solar masses as the source of energy for 7-ray 
bursts. Candidates satisfying these requirements include 
exploding yery massive stars, also terme d hypernovae or 
collapsars ( Woosley 1993 , iPaczyiiski 1998 ]) , mergers of two 
neutron stars ( Eichlcr ct al. 198S) and the merger of a 
neutron star and a black hole ( Mochkovitch ct al. 1993|) . 
Mergers may be inconsistent w ith the small offsets b e tween 
7-ray bursts and their hosts ( Bloom ct al. 1999b| , Bulik 



supernovae would often go unnoticed due to the much 
brighter 7-ray burst afterglow; therefore, we do not know 
whether most 7-ray bursts are associated with supernovae 
or just some of them. Also, the fact that GRB 980425 was 
very nearby and subluminous in gamma rays by a factor 
10 4 makes it hard to extrapolate its supernova association 
to normal 7-ray bursts without additional considerations. 
To elucidate the 7-ray burst-supernova association fur- 
ther, we examine optical and X-ray extinctions of 7-ray 
burst afterglows (Sect. ||). Then we discuss the evidence 
these extinctions provide that the majority of long 7-ray 
bursts occur in molecular clouds, with dust destruction 
explaining the unusually low extinctions (Sect. ||). We 
thus infer that neutron star and black-hole mergers are no 
longer plausible 7-ray burst progenitors, and briefly men- 
tion some further implications of our results (Sect. |j). 

2. EXTINCTION OF 7-RAY BURST AFTERGLOWS 



2.1. Optical extinction 

A paradox in associating most 7-ray bursts with explod- 
ing very massive stars is that one expects the majority of 
these to lie amidst highly absorbing molecular clouds. The 
lack of high rcstframe visual extinction has therefore led to 
considerable skepticism about the 7-ray burst-supernova 
connection. Here we reinvestigate the extinction for 7-ray 
bursts in a systematic way, by fitting an extinction model 
to all afterglows for which the required X-ray and optical 
data exist. The model function is 

F v = F vo [u/u v }- exp[-A Vl (l +z)v/v v ], (1) 

where vy and Fyo are the observer-frame V band central 
frequency and extinction-corrected V band flux, respec- 
tively, and Ay r is the rest frame visual extinction; the ex- 
tinction term is applied only to optical and infrared data. 
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et al. 1999 ). Further evidence for the relation between 
7-ray bursts and star formation has been provided by the 
fact that the brightness distribution of 7-ray bursts agrees 
well with models in which the 7-ray burst rate tracks the 
star for mation rate oyer the past 15 billion years of cosmic 



histo ry ( Totani 1997 , Wijers et al. 199S , Kommers et al 



2000). The most direct evidence relating 7-ray bursts to 



a specific type of progenitor has been the discove ry of su- 
pernova 1998b w in the error box of GRB 980425 flGalamg] 
et al. 1998b| ) and the detection of supernova-lik e light 
curves under neath two afterglow s, GRB 980326 ( Bloom 
et al. 1999a| ) and GRB 970228 flReichart 1999J , palaim 



The results of the fits are given in table H. Note the low 
optical extinctions. (Some fitted Ay,- are negative, as ex- 
pected from a fit procedure if the values are less than the 
fit errors; we did not force Ay > 0.) We have used the 
simplest possible extinction law here, A v ex v. In metal- 
rich environments, the extinction can have features, such 
as the 2200A 'bump'. We see no significant detection of 
this bump in any afterglow, despite its easy observability 
at redshifts 1-2, and therefore neglect it. 

The prompt and afterglow X rays are not significantly 
attenuated, and there is only one 7-ray burst in which 
the search for an optical afterglow resulted in upper limits 
well below the expected optical flux (GRB 970828; Groot 
et al. fl99q 4 .) Most optical non-detections are adequately 



explained by adverse observing conditions, and consistent 
with the presence of a normal afterglow. Hence there is no 
evidence from non-detections for a significant population 
of highly extincted 7-ray bursts, or for significant skewing 
of the extinction distribution due to the selection effect 
that 7-ray bursts with higher extinction are simply not 
found. 

In the above analysis, we have assumed that the intrinsic 
spectrum of the afterglow is a pure power law. The the- 
ory of afterglows allows a so-called 'cooling break' in the 
spectrum between X rays and optical, where the power- 
law slope steepens by 1/2. A comparison of the mean 
optical-to- X-ray slope, Pox, and the X-ray spectral slope, 
(3x, shows that they are indeed different, indicating that 
such a break may be present: < fix— Pox >=0.31 ± 0.05. 
Therefore we redo our fits with an extra break. The fit 
function is 

F v = F V0 [v/v v ]-W +1 M exp[-A Vr (l + z)vjv v \ 

[i/(l + z)/v c ]~ n / 2 + 1 x j " 



Ml + z)M]-"/2 + l 



(2) 



The index n controls the sharpness of the break. We use 
n — 2, which provides a smooth transition over 2 dex in 
frequency, and find little change in the fits if we vary the 
width in the range 1-3 dex. The extra datum needed to 
constrain this fit is fix- To treat this uniformly with the 
other data, we add an artificial data point with ten times 
greater frequency than the first X-ray point, and give it a 
value and error so that a fit to the two X-ray points gives 
the correct values and errors of the X-ray flux and spectral 
slope. The results of the new fits are also given in table [l| 
showing hardly any change in the extinctions. 

2.2. X-ray absorption 

Our results confirm in a more rigorous way the casual 
impression of low optical extinction in 7-ray bursts and 
therefore seem to contradict the notion of exploding mas- 
sive stars being the progenitors of 7-ray bursts. However, 
the hydrogen column densities towards the 7-ray bursts re- 
veal startlingly contrasting evidence. In table III we collect 
from the literature values of the hydrogen column density 
towards 7-ray bursts as derived from soft X-ray absorp- 
tion. To obtain these values, we subtracted the Galactic 



tx oc E~ 2 - & ( Morrison fc McCammon 1983| ), which im- 
plies that we have to multiply the foreground-subtracted 
column by (1 + z) 26 for a source at redshift z, to get the 
true column density in the rest frame of the 7-ray burst. 
(Note that while in principle the photoelectric absorption 
has richer structure due to absorption edges of individual 
atomic species, the presently available spectral resolution 
and S/N does not allow this to be discerned.) This is a 
large correction, ranging from 4 to 50 among the sample, 
and the resulting restframe hydrogen column densities are 
in the range 10 215 — 10 233 cm -2 . To emphasize the con- 
trast between the optical and X-ray results, we show in 
Figure [D the correlation between optical and X-ray extinc- 
tion. The solid curve indicates the relation between Ay 



foreground column ( Dickey fc Lockman 199C ) from the to- 
tal measured value. Then we accounted for the fact that 



and N H for the Milky Way flPrcdchl fc Schmitt 1995| ), 
so we see that the observed visual extinctions are 10-100 
times smaller than expected for the observed X-ray ab- 
sorption. 

Both types of extinction are due to heavy elements, so 
mctallicity differences cannot change the ratio. However, 
the optical/UV extinction is due to dust grains, whereas 
the X-ray extinction is due to K- and L-shell electrons of 
intermediate-mass elements (mostly C and O) and there- 
fore does not depend on whether the atoms are in a gas or 
a solid. The X-ray extinction is therefore a better measure 
of the total column density. However, converting the X-ray 
absorption to a hydrogen column density, as is customary, 
does depend on metallicity. Since the high-redshift regions 
we are probing may have lower metallicities, the true col- 
umn densities can only be larger than those we have de- 
rived. (Since the regions are very actively star forming, 
their metallicity may not be much less than in the Milky 
Way, though.) 

3. IMPLICATIONS FOR 7-RAY BURST ENVIRONMENTS 

The values of the X-ray column densities are very high, 
and typical of the c olumn densities throu gh giant molec- 
ular clouds (Fig. p| [Solomon ct al. 1987 ). This strongly 
suggests that most 7-ray bursts are located in molecular 
clouds. We now attempt to constrain the size of the clouds, 
and thereby their mass. 

First, the surprisingly low extinctions may be explained 
by recently proposed dust destruction: dust is sensitive to 
the UV and X-ray radiation from the 7-ray burst and its 
afterglow. UV and X-ray light heats the grains, and out 
to about 20 pc can evapo rate them (Waxman & Draine 
2000 , Fruchter et al. 200C ) . Waxman and Draine find that 
a prompt flash like that seen only in GRB 990123 is needed 
to muster enough UV light; however, Fruchter et al. show 
that the X-ray flux is at least as efficient in heating the 
grains. Since these come mostly from the prompt burst 
emission, the uncertain UV flashes are not needed in their 
model. They find that average 7-ray bursts can evaporate 
all dust out to 20 pc. Beyond this, dust may be shattered 
by strong grain charging. Since the effect of this on dust 
extinction properties is unclear, we shall not consider it 
here. From our low optical extinctions, we conclude that 
dust destruction must be taking place. This limits the 
bulk of the cloud to lie within about 20 pc of the 7-ray 
burst. 



follows from the fact that the afterglow radiu s after a day 
is about 0.1 pc (e.g., Wijers & Galama 1999). This firmly 
excludes any remains of the exploded star as the source 
of X-ray or optical absorption. The size of the absorbing 
cloud can also be bounded from below using absorption 
lines of Mgi in the spectra of many optical tr ansients (e.g. 



Vreeswijk et al. 2000] , [Metzger et al. 1997J ). Afterglows 



tend to have strong Mgi lines, especially relative to Mgll, 
indicating they originate in denser regions than the nor- 
mal diffuse ISM. We therefore suppose they originate in 
the same region that causes the large X-ray columns. The 
fact that this Mgi is still visible after a day means it has 
not all been ionized away. Mgi has an ionization energy 
of 7.84 eV, and photons above 13.6 eV are stopped by H 
very near the 7-ray burst. Averaging the ionization cross 
section over this energy range, weighted by the afterglow 
fluence spectrum, we find an optical depth to Mgi ionizing 
photons of 0.7Ah22.5 (assuming the solar value of [Mg/H]). 
This means that for most of the column density range, the 
optically thin limit is adequate for judging the survival of 
Mgi at the edge of the cloud. Recombination times are 
much too long to play a role. 

Integrating the 7-ray burst flux over time and over the 
same energy range, weighted by the energy-dependent 
cross section, we find that an average Mg atom would 



interc ept 2.1 x 10 4 K^V/ 2 e| /6 oe^^d/lpc)" 2 ionizing 



photons (for a typical afterglow with (3 = 0.75, where 
€b = 0.01es,_2 and e G = 0.1e e ,_i are the equiparti- 
tion fractions in the notation of Wijers & Galama 1999). 
Therefore, some of the Mg must be many parsecs from 
the 7-ray burst in order to survive. The Mg lines are 
usually rather saturated, so we can only get a lower limit 
to the Mgi column de nsity. For the case of GRB 990510 
( Vreeswijk et al. 20"00| ), this limit is 10 133 cm" 2 , which for 
normal Mg abundance is contained in less than 10 21 cm~ 2 
of total ./Vh, accounting for depiction. This is less than 
the total observed X-ray column, so only a fraction of Mg 
need remain neutral. To set a conservative lower limit, we 
shall tolerate 10 ionizations for the average Mg atom at 
the edge, implying survival of less than e -10 =5x 10 -5 
of Mg anywhere in the cloud. Then the lower limit to the 



cloud size becomes 45 El 



2/3 1/4 5/12 2/3 



'•2 



-B,-2 c e, 



e c _i pc This is quite 



sensitive to burst parameters, and often above the upper 
limit from dust destruction. Therefore, Mg lines need not 
occur in every burst, and their presence should be corre- 
lated with burst strength. 

Together, dust destruction and Mgi survival constrain 
the size of the cloud to be tens of parsecs. The cloud there- 
fore has a density of 500 AH22.5/-R20 cm -3 and a mass of 

4 X 10 5 A^H22.5-Rlo M © ( where ^1*22.5 = N H /W 22 - 5 Cm" 2 

and R20 = i?/20pc). These para meters are very much 
like t hose of giant molecular clouds ( polomon fc Edmund; 



1980n polomon et al. 1987 ). We therefore consider our 



findings strong evidence that almost all (long) 7-ray bursts 
are associated with giant molecular clouds, and therefore 
with star-forming regions. This, in turn, speaks in favor of 
massive stars rather than compact-object mergers as the 
progenitors of 7-ray bursts. 

Further predictions of the location of 7-ray bursts in 
molecular clouds are associated with absorption/scattering 



frared may produce a reradiation echo with a thermal spec- 
trum that peaks in r estframc infrared on a ti me scale of 
several tens of days ( Waxman fc Draine 2000J ). (ii) Scat- 
tering of the afterglow's light by dust outside the dust- 
vacated region may produce a scattering echo on tim e 
scales of tens to hundreds of days (Esin & Blandford 2000). 



This echo has a spectrum similar to that of the afterglow. 
Each echo can emit 10 41-42 erg/s. (iii) The far UV radia- 
tion will be absorbed by H2, causing a strong drop in the 
UV at 1650A and 1300 A, and fluorescence will result in 
restframc UV emission on time scales of days to months 
(Draine 2000). (iv) If burst radiation is collimated, it is 
likely that the later, softer emission is less collimated, en- 
abling us to see more afterglows than 7-ray bursts (e.g., 
Rhoads 1997). However, because dust is destroyed only 
along the collimated path of the initial hard radiation, 
such '7-ray burst-less afterglows' would not be visible in 
optical and near-IR from embedded sources. Only in far- 
IR, mm, and radio could the frequency of afterglows be 
significantly greater than that of 7-ray bursts. 

Our findings may also have some indirect bearing on 
the issue of the cosmic star formation history, in the fol- 
lowing sense: there has been much recent debate on the 
relative importance of UV and far-IR radiatio n in counting 
the s t ar formation rate at high redshift (e.g., Madau et al. 
1998, Barger et al. 1999). In both cases, one counts the 
location of massive, UV-producing stars, but in the far- 
IR case it is found/assumed that the majority of these are 
deeply shrouded in du st, concentrated in ultralu minous IR 
galaxies (ULIGs, e.g. [Banders fc Mirabel 1996 ). Since 7- 
ray burst radiation escapes fairly well even from ULIGs, 
7-ray burst locations might provide an unbiased sample of 
massive-star locations at redshifts 1-4. This means that 
a far-IR study of 7-ray burst host galaxies should help re- 
solve the issue of what type of host, ULIGs or UV-emitting 
smaller galaxies, are the dominant source of massive-star 
production at these redshifts. 



j-i„ „ ~i „. 



4. CONCLUSIONS 

We have examined a complete sample of 7-ray burst af- 
terglows, namely those with known redshift and X-ray col- 
umn density, for which optical to X-ray data allow a deter- 
mination of reddening. As a sample, they provide strong 
evidence for high X-ray column densities, without a single 
good exception. However, the individual measurements in 
a given source are very significant only for GRB 980703 
and GRB 980329. Therefore, good X-ray spectra are re- 
quired for more sources in order to confirm our findings 
and pin down the parameters of the clouds better. The 
size and origin of the absorbing matter are constrained by 
the low extinction, the blast wave size, and the survival of 
Mgi. High-resolution measurements of the Mg absorption 
lines to better determine the location and column density 
of the Mg absorber are needed to strengthen the lower 
limits on cloud size and mass. 

In short, we find high X-ray column densities and low 
optical extinctions for 7-ray burst afterglows, from which 
we infer that (i) Most 7-ray bursts are embedded in large 
molecular clouds, (ii) 7-ray bursts are therefore likely pro- 
duced by dying massive stars, and not by mergers of neu- 
tron stars and/or black holes, (iii) The low optical extinc- 



for the optical afterglow to escape even large clouds, (iv) 
7-ray burst host studies may help identify the dominant 
sources of high-redshift star formation. 



We are grateful to Shri Kulkarni, Daniel Reichart, James 
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supported by the Sherman Fairchild Foundation. 
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Fig. 1. — The hydrogen column density TVh versus the optical 
extinction at V bancLA^r for a number of 7-ray burst afterglows; 
data ar e from tahle h i The soli d curve is the Galactic Av-Nn re- 
lation ( Predehl &: Scnmitt 1995J ). The solid symbol is the weighted 
mean 01 the tour lowest points, assuming 50% for the unknown error 
of TVh for GRB 990123. The dashed line shows the average column 
density of a giant molecular cloud of 170Mq/pc 2 (Solomon et al 
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Table 1 

The results of fitting a single power law plus dust absorption (eq. hh and the results of fitting a broken power law plus dust absorption 
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-4.72±° i| 0.90±£j» O.SStVol 0.2/1 
-4.09±g;S 1.05±g;» 1.63±°:£ 2.0/1 
-4.94tH 6 6 0.66+^-0.20+^ 2.0/5 
-4.20+° : ° 8 9 0.86±2-°*-0.22±2:3 0.15/2 


-B.ii8:| o.46±S:S-o.t±» : s 0.3/1 
-4.38i»:iS o-so±g:ig o.i 6 +°; 2 6 7 2.1/3 

-4.81± ;« 0.38±8;g§ 0.32+°;^ 1.8/1 

-4.o6j;ro7 i.oi±g:SS ^tH 2.1/1 

-4.93lg : » 0.66±g-«-0.07±S:S* 18.8±g:« 1.9/5 
-4.09±„;:[g 0.60l°;2o O-Oi-o'.is 16-7J:?:? °-0/ 3 


0- 96 -0.i9 3 -4 0- 695 

0.99t .o? 6ig° 0.835 
1.03±g|S 200l 2 ^ 3.42 
18+. 2 . 36± 22 . 0.966 

°- 97 -o'.04 14 L6 ° 
1.12to'l 2 10tg 1.62 


980329 0402.99 
980519 0520.34 


-5.70+";^ 0.91± ;$ 2.9+^ 0.0/0 
-4.34+°^ 0.98l°;°« 0.8+.°;* 0.2/2 


-5.65t '22 0.85lo; 2 p 3.0lJ;g 0.4/0 
-4.31±° : ^ 0.92+.g : 2 J, 0.91°;^ 0.4/2 


1 63 +0 - 41 10+ 8 

± -°°_0.36 lu -4 
, t-9+0.70 4 +10 
i - OZ -0.57 -5 



